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I.  INTRODUCTION 


The  Test  and  Evaluation  Command  (TECOM)  funded  the  Experimental  Design 
and  Analysis  Branch,  SECAD,  BRL  to  develop  procedures  for  planning  durability  tests 
which  allow  testing  an  item  beyond  its  required  durability  mission  time  so  that  statisti¬ 
cal  inferences  concerning  durability  can  be  made  with  higher  confidence.  This  report 
gives  recommended  procedures  to  accomplish  this  goal,  summarizes  the  assumptions 
underlying  these  procedures,  outlines  the  theory  required  for  their  development  and 
gives  tables  and  sampling  plans  which  are  needed  to  implement  them. 

A  durability  test  is  conducted  to  collect  data  which  is  used  to  draw  inferences  con¬ 
cerning  the  probability  that  a  unit  which  is  to  be  evaluated  will  be  able  to  operate  for  a 
minimum  specified  time,  number  of  cycles,  etc.  without  failure  due  to  wearout.  If  one  is 
willing  to  agree  that  units  wear  out  gradually,  i.e.,  the  physical  properties  of  a  system 
gradually  change  with  system  usage  and  such  changes  adversely  affect  the  effectiveness 
of  the  system,  then  one  is  led  to  the  class  of  statistical  distributions  which  are  said  to 
have  an  increasing  failure  rate.  This  is  the  class  of  distributions  which  have  been  used 
for  the  work  described  in  this  report. 

The  theory  which  underlies  the  procedures  presented  in  this  report  is  described  in 
the  appendices  so  that  the  mathematically  inclined  reader  can  satisfy  his  academic 
curiosity.  The  main  body  of  the  report  is  concerned  with  application  of  the  procedures. 
The  tables  which  are  needed  to  use  the  procedures  and  listings  of  the  computer  pro¬ 
grams  which  were  used  to  construct  these  tables  are  given  in  the  appendices. 


II.  PREVIOUS  WORK 

A  durability  life  distribution  function  is  a  mathematical  equation  which  is  used  to 
calculate  the  probability  that  the  durability  life  of  a  unit  will  be  at  least  as  large  as  a 
value  which  is  of  interest  to  the  engineer  who  is  evaluating  it.  For  sake  of  consistency 
with  the  statistical  literature,  we  will  let  F(T)  be  the  durability  life  distribution  func¬ 
tion,  i.e.,  F(T)  is  the  probability  that  the  durability  life  is  less  than  or  equal  to  T  and 
R(T)  =  l-F(T)  be  the  probability  that  durability  life  is  greater  than  T.  Clearly,  the 
form  of  F(T)  depends  on  the  tasks  which  a  unit  is  required  to  perform.  For  example, 
the  F(T)  for  a  vehicle  which  is  required  to  ford  streams  and  climb  mountains  would  cer¬ 
tainly  differ  from  the  F(T)  for  the  same  vehicle  if  it  were  only  required  to  operate  on 
interstate  highways. 

There  are  many  equations  which  are  used  for  F(T).  Each  has  special  properties 
which  make  it  appropriate  for  use  with  a  certain  type  of  item.  The  two  most  commonly 
used  are  the  Weibull  distribution  and  a  special  form  of  this  distribution  which  is  known 
as  the  exponential  distribution.  The  Weibull  distribution  has  three  parameters  which 
can  be  adjusted  so  that  the  function  adequately  characterizes  the  durability  distribution 
of  many  different  types  of  equipment.  If  an  engineer  has  sufficient  knowledge  about  the 
form  of  the  durability  life  of  an  item  to  guide  him  to  a  reasonable  selection  of  the  form 
of  its  durability  distribution  and  sufficient  data  to  estimate  the  parameters  of  his 
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postulated  distribution,  then  he  should  use  this  information  to  determine  whether  the 
unit  satisfies  its  durability  requirement. 

Usually  little  is  known  about  the  functional  form  of  F(T).  Since  durability  tests  are 
almost  always  terminated  before  all  units  on  test  experience  a  durability  failure,  there  is 
seldom  sufficient  data  to  calculate  good  estimates  of  the  parameters  of  F(T),  even  if  its 
form  were  known.  This  has  led  to  the  development  of  analytical  techniques  which 
require  little  knowledge  about  the  form  of  F(T). 

The  distribution  free  procedure,  which  requires  minimal  knowledge  about  the  form 
of  F(T),  is  the  one  most  frequently  used  today.  For  this  procedure,  a  number  of  units, 
N,  are  selected  as  typical  of  the  production  design  configuration  of  the  item  to  be 
evaluated.  These  units  are  then  tested  under  a  specified  set  of  test  conditions  for  a 
durability  time  T*,  which  is  the  durability  life  requirement  for  the  item.  The  number  of 
units  which  survive  the  test,  r,  is  used  to  calculate  an  estimate  of  the  durability  of  the 
item  represented  by  the  sample.  Thus,  since  N  units  are  tested  and  r  of  them  survive 
the  test,  the  best  estimate  of  durability  is  given  r/N.  The  binomial  distribution  is  used 
to  construct  a  lower  confidence  limit  for  this  estimate.  If  one  is  willing  to  use  Bayesian 
procedures  and  postulate  a  uniform  prior  distribution  for  the  binomial  parameter,  he 
can  improve  this  estimate,  i.e.,  he  can  construct  a  lower  confidence  limit  for  durability 
which  is  larger  than  the  one  found  using  the  classical  or  frequentist  approach  to  the 
problem.  The  improvement  attained  via  this  Bayesian  approach  is  equivalent  to  having 
a  sample  size  of  N+l  with  r+1  survivors  using  the  frequentist  approach  (see  Appendix 
A). 


The  distribution  free  procedure  requires  a  large  value  of  N  if  reasonable  lower 
confidence  limits  for  durability  are  needed.  For  example,  if  the  appropriate  requirement 
document  specifies  that  a  durability  of  at  least  .80  relative  to  a  mission  time  of  5000 
miles  must  be  demonstrated  with  95%  confidence,  then  a  sample  size  of  14  units  must 
be  tested  for  time  T*  =  5000  miles  each,  with  no  durability  failures.  If  one  durability 
failure  is  to  be  allowed,  a  sample  size  of  22  units  must  be  tested  for  the  same  number  of 
miles.  If  the  item  under  evaluation  had  been  a  prototype  vehicle,  such  sample  sizes 
would  not  have  been  at  all  practical  —  sample  sizes  of  one  to  three  would  be  much  more 
reasonable. 

This  sample  size  problem  can  be  circumvented  to  some  degree  if  one  tests  for  a 
time  r  which  is  greater  than  T*  and  makes  reasonable  assumptions  concerning  the  form 
of  F(T).  This  approach  would  permit  one  to  make  the  same  inferences  about  the  dura¬ 
bility  of  a  system  relative  to  the  mission  time  T*  with  a  smaller  number  of  test  units. 
The  only  problem  associated  with  this  approach  is  the  determination  of  "reasonable” 
assumptions  concerning  the  form  of  F(T). 

Billings  (1967)  introduced  a  concept  which  he  called  a  moderately  distribution  free 
procedure.  If  F(T)  is  the  probability  distribution  for  durability  life,  then  its  first  deriva¬ 
tive  with  respect  to  T,  f(T),  is  called  the  probability  density  function  for  durability  life. 
Billings  showed  that  if  f(T)  is  an  increasing  function  of  T  during  the  test  period,  then 
for  any  number,  M,  which  is  greater  than  one,  the  relationship  F(MT)  >  MF(T)  holds. 
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He  used  this  relationship  to  develop  a  methodology  which  improves  on  the  distribution 
free  lower  confidence  limit  for  durability  when  the  actual  test  time  for  each  unit  is 
greater  than  T*.  Billings’  confidence  limits  can  be  improved  by  using  Bayesian  pro¬ 
cedures.  Though  this  work  is  mathematically  interesting,  it  has  not  been  used  by  test 
engineers  because  there  seems  to  be  little  practical  justification  for  the  assumption  that 
f(T)  is  increasing  in  T,  especially  if  the  actual  test  time  is  large  relative  to  T*.  For 
further  discussion  of  this,  see  Appendix  D. 

Barlow  and  Proschan  (1981)  consider  the  family  of  durability  probability  distribu¬ 
tion  functions  which  have  an  increasing  failure  rate  (IFR).  For  these  distributions  one 
need  only  assume  that  the  item  to  be  evaluated  has  the  property  that  its  remaining  life 
is  decreasing  with  usage,  i.e.,  things  wear  out.  They  prove  that  if  a  complex  system 
contains  components  which  independently  have  IFR  durability  distributions,  then  the 
system  has  an  increasing  failure  rate  average  (IFRA)  life  distribution.  An  IFRA  durabil¬ 
ity  distribution  has  the  property  that  the  average  of  the  failure  rate  over  time  is  increas¬ 
ing.  It  can  be  shown  that  an  IFR  durability  distribution  is  IFRA  but  an  IFRA  durabil¬ 
ity  distribution  is  not  necessarily  IFR.  Aside  from  such  generally  accepted  statements 
as,  ”90%  of  the  tire  trouble  occurs  in  the  last  10%  of  tire  life,”  the  mathematical  proofs 
given  by  Barlow  and  Proschan  adequately  justify  the  use  of  their  procedures  based  on 
IFRA  distributions  for  durability  testing  at  TECOM. 


HI.  THE  IFRA  FAMILY  OF  DURABILITY  DISTRIBUTIONS 

A  distribution  F(T)  is  said  to  have  IFRA  if  -  ln[l-F(T)]/T  is  increasing  in  T,  where 
T  is  greater  than  zero.  This  mathematical  statement,  which  leads  to  a  useful  procedure 
for  using  survival  data  for  equipment  when  test  time  exceeds  durability  requirement  to 
construct  improved  lower  confidence  bounds  for  durability,  has  a  practical  interpreta¬ 
tion.  Simply  stated,  a  system  which  has  IFRA  has  the  property  that,  even  though  there 
are  brief  periods  during  which  it  may  be  showing  a  slight  decrease  in  failure  rate  due  to 
something  like  work  hardening,  its  failure  rate  is  increasing  when  averaged  over  time. 
Mathematical  developments  which  follow  from  the  IFRA  assumption  are  given  in 
Appendix  B. 

The  concept  of  IFRA  evolved  from  the  fact  that  most  complex  systems  have  com¬ 
ponents  that  wear  out,  i.e.,  are  IFR.  Barlow  and  Proschan  (1967)  proved  that  any  rea¬ 
sonable  system  built  from  components  having  either  exponential  or  IFR  durability  dis¬ 
tributions  will  have  an  IFRA  distribution.  It  can  also  be  shown  that  if  each  of  the 
independent  components  of  a  system  has  an  IFRA  distribution  then  the  system  has  an 
IFRA  distribution.* 

Barlow  and  Proschan  (1981)  give  another  interesting  justification  for  the  use  of  the 
IFRA  distributions.  They  prove  that  an  IFRA  distribution  arises  naturally  when  ran¬ 
dom  shocks  occur,  each  independently  causing  random  damage  to  a  device,  the  damages 


•Theorem  2.6,  Barlow  and  Proschan  (1081). 
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accumulating  until  a  critical  threshold  is  exceeded,  at  which  time  the  device  fails.  This 
is  called  the  cumulative  damage  model. 

Thus,  the  IFRA  distributions  seem  to  be  a  natural  for  durability  testing  and  are 
totally  consistent  with  the  usual  assumption  that  a  unit  which  has  been  in  the  field  for 
several  years  is  more  likely  to  wear  out  tomorrow  than  one  which  has  been  in  the  field 
for  several  months. 

All  of  the  mathematical  models  which  underlie  the  procedures  given  in  the 
remainder  of  this  report  are  based  on  the  assumption  that  the  appropriate  durability 
probability  distributions  are  IFRA. 


IV.  DESIGN  OF  EXPERIMENTS 

The  durability  requirement  for  an  item,  given  in  terms  of  some  measure  appropri¬ 
ate  to  the  item,  e.g.,  hours,  cycles,  miles,  rounds,  etc.,  should  specify  T*,  the  required 
durability  life  for  the  item,  R*,  the  required  durability  relative  to  T*;  and  a  level  of 
confidence  at  which  R*  must  be  demonstrated.  Implicit  in  this  requirement  is  the  life 
profile  for  the  item  and  the  definition  of  a  durability  failure.  Further  discussion  of  these 
points  can  be  found  in  Test  Operations  Procedure  1-2-502  (Reference  7). 

If  a  durability  test  is  to  demonstrate  the  durability  life  of  an  item,  two  things  must 
be  determined,  the  number  of  units  to  be  tested,  N,  and  the  time  on  test  for  each  unit, 
r.  Procedures  for  selecting  N  and  r  when  the  durability  life  of  an  item  is  IFRA  are 
given  below.  A  brief  discussion  of  the  theory  underlying  the  procedure  is  given  in 
Appendix  B,  and  a  more  comprehensive  derivation  is  given  in  Barlow  and  Gupta  (1967). 

Given  a,  the  confidence  coefficient,  T*,  and  R*,  Table  F-l  can  be  used  to  select  N, 
the  number  of  units  to  put  on  test,  and  r,  the  test  time  for  each  surviving  unit,  so  that 
a  100(1  -  a)%  lower  confidence  limit  for  R(T*),  denoted  as  R(a,T*),  will  be  at  least  R* 
provided  that  the  test  time  is  r  and  r  units  survive.  It  should  be  pointed  out  that  each 
unit  must  remain  on  test  until  it  either  fails  or  has  operated  for  time  r. 

Since  r  must  be  greater  than  T*,  we  can  write  r  =  MT*,  i.e.,  r/T*  =  M,  where  M 
is  a  number  greater  than  one.  This  permits  construction  of  a  concise  table  which  has 
general  application. 

The  use  of  Table  F-l  to  select  a  sample  size  is  best  illustrated  with  an  example. 
Suppose  that  there  is  a  requirement  that  a  new  gun  tube  must  demonstrate  a  durability 
of  at  least  .80  with  95%  confidence  when  firing  2400  rounds  of  high  velocity  K.E. 
ammunition.  That  is  to  say,  the  requirement  is  to  demonstrate  with  95%  confidence 
that  the  probability  that  tube  life  exceeds  2400  rounds  is  at  least  .80.  Here  a  =  .05, 
T*  =  2400  and  R*  =  .80.  Table  F-l  gives  sample  sizes  such  that,  if  N  units  are  each 
tested  for  2400M  rounds  and  C  =  N  -  r  units  fail,  then  R(.05,  2400)  is  at  least  .80  with 
at  least  95%  confidence.  Of  course,  it  is  assumed  that  the  R(T)  for  the  gun  tube  is 
IFRA. 
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For  this  example  a,  T*  and  R(+2400)  are  all  fixed  and  the  test  engineer  is  free  to 
select  N,  M  and  C  in  such  a  way  that  the  durability  requirements  are  satisfied  and,  at 
the  same  time,  the  test  can  be  run  as  cheaply  as  possible.  If  no  durability  failures  are 
allowed,  i.e.,  C  =  0  and  M  =  1.7,  so  that  all  units  on  test  must  fire  MT*  =  1.7x2400 
=  4080  rounds  without  a  durability  failure,  then  a  sample  of  N  =  8  gun  tubes  is 
required.  This  means  that  a  total  of  8x4080  =  32640  rounds  would  be  fired.  It  is  of 
interest  to  note  that  if  the  distribution  free  method  were  to  be  used  for  this  demonstra¬ 
tion,  2400  rounds  would  be  fired  from  each  of  14  gun  tubes.  That  is,  the  number  of 
rounds  fired  would  be  14x2400  =  33600,  only  960  rounds  more  than  with  the  IFRA 
design.  But  six  additional  gun  tubes  would  be  required  and  this  would  also  increase  the 
cost  of  the  test. 

Quite  often  the  number  of  units  available  for  test  is  fixed.  If  this  were  the  case  in 
the  above  example  and  only  N  =  3  gun  tubes  were  available  for  testing,  then  Table  F-l 
reveals  that,  if  C  =  0,  i.e.,  no  durability  failures  would  be  allowed,  then  M  must  be 
greater  than  3,  which  is  not  recommended  because  extrapolation  back  from  more  than 
three  times  the  durability  requirement  is  mathematically  undesirable.  In  this  case,  the 
test  engineer  would  most  likely  request  that  the  requirement  for  durability  demonstra¬ 
tion  with  95%  confidence  be  changed  to  durability  demonstration  with  80%  confidence. 
The  value  of  M  from  Table  F-l  for  this  requirement  is  2.6.  Then  2.6x2400  =  6240 
rounds  would  be  fired  from  each  of  the  three  gun  tubes,  and  if  there  were  no  durability 
failures,  then  the  test  would  have  demonstrated  that  the  durability  for  a  2100  round 
mission  is  at  least  .80  with  82%  confidence. 


V.  ESTIMATION  OF  DISABILITY 

The  previous  section  gave  procedures  for  selecting  sample  sizes  so  that,  if  C  or  less 
failures  were  observed  in  a  sample  of  N  test  units,  a  100(l-a)%  lower  confidence  limit 
for  durability  would  be  at  least  as  large  as  a  stated  requirement.  After  the  test  has  been 
completed  and  the  resultant  data  becomes  available  for  analysis,  a  better,  i.e.,  higher, 
lower  confidence  limit  for  durability  can  usually  be  constructed.  Also,  point  estimates  of 
durability  for  any  stated  time  requirement  and  time  for  which  a  specified  durability  is 
predicted  can  be  generated  when  C  >  0.  As  before,  the  IFRA  durability  distribution 
-assumption  is  made.  Methodology  for  doing  the  above  is  given  in  this  section. 

When  a  durability  test  has  been  completed  the  evaluator  knows  the  number  of 
units  which  were  on  test,  the  number  of  units  which  failed,  the  time  at  which  each 
failure  occurred  and  the  time  at  which  the  experiment  was  terminated.  If  the  termina¬ 
tion  time  is  the  same  as  the  required  durability  mission,  time,  i.e.,  T*  =  r,  then  the  best 
point  estimate  of  R(T*)  is  given  by  R(T*)  =  r/N,  no  matter  what  form  R(T)  may  have. 

When  t  is  greater  than  T*  and  R(T)  is  EFRA,  we  know  that  R(T)  is  bounded  below 
by  the  exponential  durability  distribution.  That  is  to  say,  when  0  <  T*  <  r  then 
R(T)  >  e~T*/^,  where  0  is  an  appropriate  estimate  of  the  mean  of  the  bounding  distribu¬ 
tion.  The  maximum  likelihood  estimate  for  6$  is 


jt 

"  See  Appendix  C  for  the  derivation  of  this  result. 
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> 


(1) 


0  = 


C 


E  Ts  +  (N  -  C )T 

i=l 


C 

where  ^  T;  is  the  total  time  on  test  for  units  which  fail  and  (N  -  C)r  is  the  total  time 
i=l 

on  test  for  units  which  survive.  Thus  0  is  the  total  test  time  for  all  units  divided  by  the 
number  of  durability  failures. 

It  should  be  pointed  out  that  6,  like  many  maximum  likelihood  estimates,  is  a 
biased  estimate  of  0.  How'ever,  in  this  case  the  small  positive  bias  seems  desirable 
because  it  tends  to  compensate  for  the  conservative  bound  for  R(T*)  given  by  the 
exponential  durability  distribution. 

For  example,  suppose  four  units  were  on  test,  one  failed  at  1095  hours  and  the 
remaining  three  were  still  operating  at  the  test  termination  time  of  2300  hours.  Here 

A  1  C 

N=4,  C=l,  t  =  2300  and  9  =  -A-  [  £  T;  +  (N  -  C)r  ]  =  1095  +  (3)2300  =  7995 

C  i=1 

hours.  If  T*  =  1800  hours  then  we  would  infer  that  R(T*)  =  R(1800)  =  e  180%9!)5  = 

.80.  However,  if  instead  we  have  T*  =  1150  hours,  then  we  would  infer  that  R(1150)  — 
e-l  150/7995  _  g7 


Durability  test  data  can  also  be  used  to  construct  lower  confidence  limits  for  dura¬ 
bility.  This  is  accomplished  by  constructing  a  lower  confidence  limit  for  the  exponential 
parameter,  9,  and  using  the  exponential  bound  for  the  IFRA  class  of  distributions.  This 
requires  knowledge  of  the  sampling  distribution  of  9. 

The  exact  distribution  of  9  was  derived  by  Bartholomew  (19G3)  for  the  case  in 
which  C  >  0.  Unfortunately,  the  function  derived  by  Bartholomew  is  mathematically 
complex  and  can  not  be  expressed  in  closed  form.  A  number  of  approximations  for  the 
distribution  of  0  have  been  derived  to  circumvent  this  problem.  If  the  test  is  terminated 
at  a  fixed  time,  then  we  suggest  using  the  approximation  that  2C0/9  has  the  chi-square 
distribution  with  2C+1#  degrees  of  freedom.  This  leads  to  the  approximate 
100(1  -  a)%  lower  confidence  limit  for  0 : 

0L  =  2C0/X\.Q  (2C+1)  .  (2) 


If  C=0,  we  use  the  binomial  distribution  to  construct  a  lower  confidence  limit  for 
durability  at  time  r,  set  this  equal  to  R(a,r),  solve  for  an(l  use  this  solution  to  find 
R(a,T*).  A  lower  100(l-a)%  confidence  limit  for  the  binomial  parameter  P  is  given  by 
the  solution  of  the  equation 

E  (x  )  Pl*  (1-Pl)N~x  =  a  (•I) 

x=r 


"  There  are  many  approximations  for  this  confidence  limit.  The  one  which  is  most  common  involves  the  U3e  of  a  cbi-square 
variate  with  2C+2  degrees  of  freedom  Preliminary  results  from  our  study  of  the  problem  suggest  that  a  chi-square  with  2C+1  d  f  , 
which  has  been  considered  by  others,  may  be  the  best  chi-square  approximation.  At  this  time  we  recommend  using  the  chi-square 
approximation  with  2C+1  d.f.  because  it  gave  confidence  limits  which  are  closest  to  the  exact  limits  in  a  recent  monte-carlo  study. 
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for  P l-  In  this  equation  a  is  the  confidence  coefficient,  r  is  the  number  of  units  surviv¬ 
ing  when  the  durability  test  is  terminated,  N  is  the  number  of  units  on  test  and  PL  is 
the  lower  confidence  for  durability,  r  =  N  when  C  =  0  so  that  the  equation  becomes 

PlN  =  «  (4) 

which  has  the  solution 

PL  =  <»1/N  ■  (5) 

But,  when  we  use  the  exponential  bound  for  the  IFRA  class  of  distributions 
Pl  =  R(o,t)  =  e~r^L  =  a!/N.  This  equation  has  the  solution 

&L  =  -Nr/ln(or) .  (6) 

Substituting  this  solution  in  our  equation  for  R(a,T*)  leads  to 

T«In(or)  T* 

R(a,T*)  =  e  Nr  =  e1"^)  (7) 

T«  1 

__  a  Nr  =a  NM 

•f 

where  M  =  —  is  the  ratio  of  test  termination  time  to  durability  mission  requirement 

time.  Of  course,  our  requirement  that  r  must  be  greater  than  or  equal  to  T*  dictates 
that  M  must  be  greater  than  one.  The  above  turns  out  to  be  a  reasonably  good  approx¬ 
imation  which  is  adequate  for  a  back-of~the-envelope  analysb  in  many  cases. 

Fortunately,  the  computing  power  required  to  numerically  integrate  the  exact  pro¬ 
bability  density  function  of  9  is  now  generally  available  and  exact  lower  confidence  lim¬ 
its  for  9  can  be  constructed.  A  computer  program  for  this  purpose  has  been  written  and 
is  available  on  the  BRL  CYBER  computer  for  general  use.  A  listing  of  this  program  is 
given  in  Appendix  E. 

Given  0^,  the  lower  confidence  limit  for  9,  then  a  bound  for  the  lower  confidence 
limit  for  durability  for  time  T*  is  given  by 

R(a,T»)  =  e-T-/'1 .  (8) 


In  the  example  above  we  had  r  =  2300,  N  =  4,  C  =  1  and  0  =  7995.  The  exact 
90%  lower  confidence  limit  for  9  is  0^  —  3053  so  that  a  90%  lower  confidence  limit  for 
durability  for  time  T*  is 


R(.10,T*)  =  e~T*/3053  .  (9) 

For  sake  of  comparison,  we  calculate  the  approximate  lower  confidence  limit  for  9  to  be 

2C0  2(1)(7QQ5) 


&L  = 


X  .oo  (2C+1) 


6.25 


2558. 


(10) 
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These  exact  and  approximate  bounds,  R(o;,T*),  were  calculated  for  several  values  of 
T*  and  are  shown  in  Table  1.  Point  estimates  are  included  in  the  table  as  reference 
points. 

TABLE  1. 

Point  and  interval  estimates  of  durability;  four 
units  on  test,  one  failed  at  1095  hours.  Test 
terminated  at  2300  with  three  units  surviving. 


T* 

R(T*) 

R(.10,T*) 

R(.10,T*) 

ERROR 

500 

.94 

.85 

.82 

-.03 

1000 

.88 

.72 

.68 

-.04 

1500 

.83 

.61 

.56 

-.05 

2000 

.78 

.52 

.46 

-.06 

2500 

.73 

.44 

.38 

-.06 

Inspection  of  Table  1  reveals  that,  for  this  example,  the  chi-square  approximation 
compares  favorably,  i.e.,  the  error  relative  to  the  exact  value  is  small.  This  table  also 
shows  that  the  approximation  is  conservative,  i.e.,  R(.10,T*)  <  R(.10,T*).  These  are 
not  general  properties  of  the  approximation.  They  apply  only  to  this  specific  example 
and  no  attempt  should  be  made  to  extend  them  to  other  problems.  We  are  now  con¬ 
ducting  a  study  which  should  identify  the  best  procedure  for  constructing  approximate 
lower  confidence  limits  for  the  exponential  parameter  when  durability  test  data  is  avail¬ 
able.  We  are  also  exploring  the  possibility  of  constructing  charts  and/or  tables  which 
can  be  used  to  find  exact  lower  confidence  limits  for  the  exponential  parameter  without 
recourse  to  a  large  computer. 

Both  the  exact  and  approximate  confidence  given  above  exist  only  if  one  or  more 
durability  failures  occur  during  a  test.  This  is  a  result  of  the  fact  that  0  is  not  defined 
f°r  C=0.  If  C=0  then  a  point  estimate  for  durability  exists  only  if  T*  =  r  and,  in  this 
case,  it  is  unrealistically  one. 

If,  as  in  the  previous  example,  four  units  were  on  test  and  the  test  was  terminated 
after  2300  hours  with  all  units  surviving,  then,  for  T*  =  1150  we  have  M  =  2  and 

j_ 

R(.10,1150)  =  (.l)8  =  .7499.  (11) 
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VI.  TABLE  USAGE  WITH  EXAMPLES 

The  tables  given  in  Appendices  F,  G  and  H  are  appropriate  for  items  whose  dura¬ 
bility  probability  distribution  functions  have  increasing  failure  rate  average  (IFRA).  An 
IFRA  durability  distribution  has  the  property  that  the  average  of  the  failure  rate  over 
time  is  increasing.  The  IFRA  distributions  seem  to  be  appropriate  for  durability  testing 
since  it  seems  natural  that  the  likelihood  of  a  failure  increases  as  the  amount  of  testing 
increases  and  the  items  wear  out.  (See  Section  III  for  more  details.) 

A.  Table  F-l. 

Table  F-l  can  be  used  to  determine  the  number  of  items  needed  for  a  durability 
test  when  the  following  items  are  specified:  durability  requirement,  number  of  failures 
allowed,  desired  confidence  level  and  M.  M  is  defined  as:  M  =  test  time/durability  cri¬ 
terion.  Time  is  a  generic  term  used  to  represent  miles,  cycles,  number  of  rounds,  actual 
time,  or  any  other  measure  of  operation  sequence. 

The  tables  have  been  generated  for  the  following  values: 


Confidence  Level: 
Durability: 

Number  of  Failures: 
M: 


.95,  .90,  .80,  .70,  .60,  .50 
.50,  .60,  .70,  .80,  .90,  .95 
0,  1,  2,  3,  4,  5,  10,  25 
1.0  to  2.0  by  0.1 
2.0  to  3.0  by  0.2  . 


Suppose  that  we  are  interested  in  demonstrating  that  a  truck  has  a  durability  pro¬ 
bability  of  0.80  after  being  driven  for  5000  miles  at  a  confidence  level  of  0.90.  We 
decide  that  we  will  accept  one  failure,  given  that  the  other  trucks  have  been  driven  8000 
miles  each.  How  many  trucks  do  we  need  to  test? 

For  this  example: 

_ test  time  _ r_  _  8000  __  j  g 

durability  criterion  T*  5000 

R*  =  durability  =  0.80 

Confidence  Level  =  1  -  a  =  0.90 

Number  of  Failures  (C)  =  1  . 

Go  to  the  10th  page  of  Table  F-l.  Under  the  Durability  column  heading  of  .80,  go 
to  the  row  for  M  =  1.6  and  read  the  values  12  and  .92.  This  indicates  that  if  12  trucks 
are  put  on  test  for  8000  miles  each  and  if  at  most  one  failure  occurs,  then  our  durability 
is  0.80  with  confidence  of  at  least  0.92.  Due  to  the  discreteness  of  the  data,  the  true 
confidence  level  is  higher  than  the  stated  confidence  level. 
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Let’s  assume  that  it  is  learned  that  we  cannot  get  12  trucks  but  the  durability 
requirement  has  been  reduced  to  0.70.  Then  we  can  move  to  the  left  in  the  same  row 
and  read  under  .70,  8  and  .93.  This  indicates  that  if  8  trucks  are  put  on  test  for  8000 
miles  each  and  no  more  than  one  failure  occurs,  then  we  can  state  that  our  durability  is 
0.70  with  confidence  of  at  least  0.93. 

B.  Table  G-l. 


units^fe^  ^  ^  ^  US6<*  t0  determine  t*le  normalized  test  time  for  surviving 


M  = 


test  time 


durability  criterion 


r 

T*  . 


(12) 


In  order  to  find  M,  the  following  items  need  to  be  specified:  confidence  level,  durability 
number  of  failures  allowed  and  number  of  test  units.  ’ 


The  tables  have  been  generated  for  the  following  values: 


Confidence  Level: 
Durability: 

Number  of  Failures: 
Number  Test  Items: 


.50,  .60,  .70,  .75,  .80,  .85,  .90,  .95,  .99 
.50,  .60,  .70,  .75,  .80,  .85,  .90,  .95 

0,  1 

2  thru  10. 


Suppose  we  have  four  tanks  available  for  testing  the  durability  of  the  gun  tube. 

inftnSPeClfijCatlODi.S  f°r  the  gUn  tube  state  that  the  durability  should  be  .70  after  firing 
iOOO  rounds  with  no  more  than  one  failure  with  confidence  of  .80.  How  many  rounds 

should  we  plan  on  shooting  from  each  tank  gun  tube  in  order  to  meet  the  requirements? 
For  this  example, 

N  ==  number  of  units  on  test  =  4 
C  =  number  of  failures  allowed  =  1 
R*  =  durability  =  0.70 
Confidence  Level  =  1  -  a  =  0.80 


M  = 


test  time 


durability  criterion 


T 

tT 


1000 


(13) 


Go  to  the  12th  page  of  Table  G-l.  Under  the  Durability  Column  heading  of  0.70 
go  to  the  row  for  CONF  =  .80  and  read  M  =  2.449.  Then 


t  =  1000  •  M  =  2449 


(14) 


This  indicates  that  2449  rounds  need  to  be  fired  from  each  gun  tube,  with  no  more 
than  one  failure,  in  order  to  meet  the  specified  conditions. 
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C.  Table  H-l 


Table  H-l  can  be  used  to  determine  the  sample  size  required  to  demonstrate  a 
mean  durability  (i.e.,  the  required  durability  life)  with  a  certain  minimum  confidence 
level.  In  order  to  use  the  tables  the  following  items  need  to  be  specified.  The 
confidence  level,  the  number  of  allowable  failures  and  two  of  the  following:  M,  r  and  T*. 


Thus,  if  any  two  are  known,  the  third  can  be  calculated. 


(15) 


The  tables  have  been  generated  for  the  following  values: 


Confidence  Level: 
Number  of  Failures: 
M: 


.95,  .90,  .80,  .70,  .60,  .50 
0  thru  10 
1.1  to  2.0  by  0.1 
2.0  to  3.0  by  0.2  . 


Suppose  that  we  wish  to  demonstrate  that  the  mean  durability  of  combat  boots  is 
10  months  with  95%  confidence.  We  plan  on  testing  the  boots  for  12  months  and  are 
willing  to  allow  up  to  six  failures. 

For  this  example: 


n .  test  time  12  ,  _ 

M  =  — - — - : — : —  =  —  =  1.2 

durability  criterion  10 

Confidence  Level  =  1  -  a  =  0.95 
C  =  number  of  allowable  failures  =  6  . 

Go  to  the  second  page  of  Table  H-l.  Under  the  C=6  column,  go  to  the  row  for 
M  =  1.2  and  read  N  =  35  and  true  confidence  =  .95.  (Due  to  the  discreteness  of  the 
data,  the  true  confidence  may  be  higher  than  the  stated  confidence  level.)  This  indi¬ 
cates  that  if  we  put  35  pairs  of  combat  boots  on  test  for  12  months  and  get  no  more 
than  6  failures,  the  mean  durability  is  equal  to  or  greater  than  10  months  with  at  least  a 
95%  confidence  level. 

Instead  of  determining  N,  Table  H-l  could  be  used  to  determine  test  time  if  the  fol¬ 
lowing  items  were  specified:  C,  1  -  a,  durability  criterion  and  N.  For  example,  let’s  set 
up  another  boot  test.  Say  N  =  20,  C  =  8,  1  -  a  =  .90  and  we  want  the  mean  durabil¬ 
ity  equal  to  or  greater  than  12  months.  Go  to  the  fourth  page  of  Table  H-l.  Under  the 
C=8  column,  go  down  the  column  until  N=20  is  found,  then  read  the  M  associated 
with  this  value  M=1.5.  We  can  then  calculate  test  time  =  M  •  durability  criterion  == 
1.5  (12)  =  18. 
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Thus,  if  we  put  20  pairs  of  boots  on  test  for  18  months  and  observe  no  more  than 
eight  failures,  then  the  mean  durability  is  equal  to  or  greater  than  12  with  at  least  a 
92%  confidence  level. 


VH.  ACKNOWLEDGMENT 

The  authors  wish  to  express  their  appreciation  to  William  E.  Baker,  BRL  and  Larry 
L.  West,  TECOM  for  their  helpful  comments  provided  after  reviewing  the  report.  Our 
thanks  also  go  to  Annette  L.  Wiseman  for  the  many  typed  drafts  and  for  the  final  typ¬ 
ing  of  the  report. 


20 


REFERENCES 


1.  Billings,  Michael  G.,  ”A  Moderately  Distribution  Free  Approach  to  Reliability  Esti¬ 
mation  Based  on  the  First  Order  Statistic,”  Proceedings  of  the  Twelfth  Conference 
on  the  Design  of  Experiments  in  Army  Research  Development  and  Testing,  ARO-D 
Report  67-2. 

2.  Barlow,  Richard  E.,  Proschan , Frank,  Statistical  Theory  of  Reliability  and  Life 
Testing.  Silver  Spring,  MD  ,  1975. 

3.  Barlow,  Richard  E.,  Proschan, Frank,  "Tolerance  and  Confidence  Limits  for  Classes 
of  Distributions  Based  on  Failure  Rate,”  Annals  of  Mathematical  Statistics.  Vol.  37, 
No.  6,  pp.  1593-1601,  December  1966. 

4.  Barlow,  Richard  E.,  Proschan,  Frank,  "Inequalities  for  Linear  Combinations  of 
Order  Statistics  from  Restricted  F amilies,”  Annals  of  Mathematical  Statistics.  Vol. 
37,  No.  6,  pp.  1574-1592,  December  1966. 

5.  Billings,  Michael  G.,  "Improving  Binomial  Reliability  Estimates  -  A  Moderately 
Distribution  Free  Technique  for  Small  Samples  Reliability  Estimation,”  Proceedings 
of  the  Thirteenth  Conference  on  the  Design  of  Experiments  in  Army  Research 
Development  and  Testing,  ARO-D  Report  68-2. 

6.  Barlow,  Richard  E.,  Gupta,  Shanti  S.,  "Distribution  -  Free  Life  Test  Sampling 
Plans,”  Technometrics.  Vol.  8,  No.  4,  November  1966. 

7.  US  Army  Test  and  Evaluation  Command  Engineering  Test  -  Common  Test 
Operation  Procedures,  AMSTE-RP-702-100,  Test  Operations  Procedure  1-2-502, 
AD  759770,  14  September  1972. 

8.  Hagan,  John,  "An  Investigation  of  the  Applicability  of  Some  Moderately 
Distribution-Free  Methods  to  Durability  Tests,”  USATECOM  No.  9-00-001-000- 
020,  5  November  1970. 

9.  "Test  Design  Plan  for  DT  II  of  the  MK19  MOD3  40mm  (Grenade)  Machine  Gun,” 
US  Army  Material  Systems  Analysis  Activity,  Aberdeen  Proving  Ground,  MD, 
January  1983. 

10.  Bartholomew,  D.J.,  "The  Sampling  Distribution  of  an  Estimation  Arising  in  Life 
Testing,”  Technometrics.  Vol.  5,  No.  3,  pp.  361-374,  August  1963. 


21 


APPENDIX  A.  DISTRIBUTION  FREE  (BINOMIAL)  CONFIDENCE  LIMITS 


23 


APPENDIX  A.  DISTRIBUTION  FREE  (BINOMIAL)  CONFIDENCE  LIMITS 


In  a  test  to  determine  the  durability  of  an  item,  N  units  are  tested  for  a  required 
durability  life,  T*,  and  the  number  of  units  which  survive,  r,  is  noted.  The  probability 
that  the  durability  life  of  a  unit  is  less  than  or  equal  to  T  is  given  by  the  function  F(T). 
Thus,  the  probability  that  the  durability  life  satisfies  the  requirement  is 
1  -  F(T*)  =  R(T*).  If  no  assumption  concerning  the  mathematical  form  of  R(T)  is 
made,  a  lower  confidence  limit  for  durability,  R(T*),  of  a  unit  can  be  constructed  using 
the  binomial  distribution.  Such  a  confidence  interval  is  called  "distribution  free” 
because  it  does  not  depend  on  the  distribution  of  T.  That  is,  it  is  a  valid  interval  no 
matter  what  the  function  R(T)  may  be. 

Two  methods  are  commonly  used  to  construct  a  distribution  free  lower  confidence 
limit  for  R(T).  These  are  the  classical  or  frequentist  approach  and  the  Bayesian 
approach.  Both  methods  are  based  on  the  probability  statement, 

P  {R(T*)  >  R(a,T*)  }  =  1  -  a  ,  (Al) 

where  the  confidence  coefficient,  a,  satisfies  the  condition  0  <  a  <  .5.  Thus,  if  a  =  .05 
the  probability  that  the  durability  of  a  unit  is  greater  than  or  equal  to  R(.05,T*)  is  .95. 

For  the  frequentist  approach  R(T*)  is  considered  to  be  an  unknown  constant  so 
that  the  probability  that  exactly  r  units  in  a  sample  of  N  have  durability  life  greater 
than  or  equal  to  T*  is  given  by  the  binomial  expression 

(^)  R(T*)r  (1  -  R(T*))N~r  (A2) 

and  the  probability  that  r  or  more  units  have  durability  life  greater  than  or  equal  to  T* 
is  given  by 

E  (1jI)  R(T*)‘ (1  -  R(T*))N~' .  (A3) 

i— r 


A  100(1  -a)%  lower  confidence  limit  for  R(T*)  is  found  by  solving  the  equation 

E  (?)  RKT.ji  (1  -R(a,T.))Ni  =  a  (A4) 

i=r 

for  R(a,T*).  This  leads  to  the  probability  statement,  P{R(T*)  <  R(a,T*)}  =o.  Here 
we  claim  that  if  R(T*)  is  less  than  or  equal  to  R(a,T*),  it  would  be  most  unusual  for  r 
or  more  units  to  survive  the  test.  For  example,  if  a  =  .05,  the  probability  that  r  or 
more  units  would  survive  the  test  if  R(T*)  <  R(a,T*)  is  less  than  or  equal  to  .05. 

For  the  Bayesian  approach  to  the  problem  R(T)  is  assumed  to  be  a  random  variable 
and  a  prior  distribution  is  specified  for  it.  The  beta  distribution  is  usually  chosen  for 
the  Bayesian  prior  distribution  of  the  binomial  parameter  because  of  its  desirable 
mathematical  properties.  Like  the  binomial  parameter,  the  beta  random  variable  can 
take  on  only  values  between  zero  and  one;  the  distribution  has  two  parameters,  A  and 
B,  which  can  be  chosen  so  that  it  can  take  on  a  wide  spectrum  of  shapes;  and  the 
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mathematics  required  for  constructing  confidence  limits  are  relatively  simple. 


In  the  absence  of  an  understanding  of  the  prior  distribution  of  R(T*),  the  uniform 
distribution  is  usually  assumed.  This  is  a  special  form  of  the  beta  distribution  with 
parameters  A  =  B  =  1.  This  implies  that  R(T*)  is  equally  likely  to  take  on  any  value 
between  zero  and  one,  i.e., 

f,  (R(T*))  =  1  ;  o  <  R(T*)  <  1 
=  0  ;  otherwise. 

Similar  to  the  frequentist  approach 

f2  (r  |  R(T*))  =  ( N)  R(T.)'  (1  -  R(T.))N-' ,  r  =  0,1 . N  (A6) 

is  the  conditional  distribution  of  r  given  R(T*).  It  follows  that  the  joint  distribution  of 
r  and  R(T0)  is  given  by 

f3  (r,R(T*))  =  h  (R(T*))  f2  (r  |  R(T*))  (A7) 

—  (?)  R(T*)r  U  “  R(T*)N*r  0  <  r(x*)  <  1 

r  =  0,  1,  ...,  N 

=  0  Otherwise  . 

The  marginal  distribution  of  r  is  found  by  integrating  equation  (A7)  with  respect  to 
R(T*),  i.e., 

MO  =  J  ( ? )  R(T*)r  ( 1  -  R(T.))n"'  dR(T.)  (A8) 


1 


N+l 


;  r  =  0,  1,  ...,  N  . 


It  is  interesting  to  note  that  we  started  with  a  continuous  uniform  marginal  distribution 
for  R(T*)  and  found  that  r  has  a  discrete  uniform  marginal  distribution. 


Bayes  theorem  specifies  that  the  conditional  distribution  of  R(T*)  given  r  is 

fS  (r.R(T*)) 


4  (R(T*)  I  r) 


Mr) 


(A9) 


(N)  R(T.)' (1  -  R(T.))N- 


N+l 


=  r  (,+i)  mnI+1)  r(t*)' (1-r(t*))n-' 


which  is  a  beta  distribution  with  parameters  A  =  r+1  and  B  =  N-r+1.  A  Bayesian 
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100(1  -  a)%  lower  confidence  limit  for  R(T*)  is  given  by  the  solution  of  the  integral 
equation 


R(o,T*) 


/ 


r  (n+2) 

r  (r  +  1)  r  (N-r+l) 


R(T*)r  (1  -  R(T*))N-r  =  a  . 


(A10) 


If  we  integrate  the  left  hand  side  of  the  equation  (A10)  by  parts,  we  find  that 
R(«,T*) 


J 


r  (N+2) 


0  r  (r  +  1)  r  (N-r+l) 
N+1  'N+l 


R(T*)r  (1  -  R(T*)) 


N-r 


(All) 


=  E  f1  i  )  Rfa.T*)1  (1  -  R(of,T*))N_i  =  Of . 

i=r+l 


Comparing  equation  (A4)  with  equation  (All)  reveals  that  the  Bayesian  lower 
confidence  limit  for  R(T*)  when  N  units  are  tested  and  r  units  survive  is  the  same  as  the 
frequentist  lower  confidence  limit  for  R(T*)  when  N+l  units  are  tested  and  r+1  units 
survive.  Thus,  the  additional  assumptions  for  the  Bayesian  analysis  with  a  uniform 
prior  distribution  for  R(T*)  are  equivalent  to  an  increased  frequentist  sample  size  of  one 
unit  which  survives. 
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APPENDIX  B.  INCREASING  FAILURE  RATE  AVERAGE  SYSTEMS 

The  failure  rate  for  a  durability  distribution  is  defined  to  be 

r,T»  =  7%)  ■  <B1> 

Clearly  r(T)  AT  is  the  probability  that  an  item  which  has  survived  to  time  T  will 
experience  a  durability  failure  at  some  time  between  T  and  T  +  AT.  For  example,  if 
f(T)  =  Xe  XT,  i.e.,  the  durability  life  has  the  exponential  distribution,  then  R(T)  =  e-XT 
and  r(T)  AT  =  X  AT.  In  this  case,  the  probability  of  a  durability  failure  in  a  time 
interval  of  length  AT  depends  only  on  AT  and  not  on  the  time  T.  Thus,  items  with  an 
exponential  durability  life  don’t  wear  out.  Their  failure  rate  is  constant  over  time. 
Such  items  tend  to  operate  without  a  durability  failure  until  it  experiences  a  random 
external  pulse  or  shock  which  induces  failure. 

A  durability  distribution  is  said  to  have  an  increasing  failure  rate  (IFR)  if 
r(T)  <  r(T  +  AT)  for  all  T  and  positive  AT.  Any  item  which  wears  out  has  IFR,  i.e., 
the  older  it  is  the  more  likely  it  is  to  experience  a  durability  failure  during  the  next 
hour,  mile,  cycle,  etc. 

A  durability  distribution  is  defined  to  have  an  increasing  failure  rate  average 
(IFRA)  if  the  average  value  of  r(T)  increases  over  time.  The  average  value  of  r(T)  over 
the  interval  from  time  0  to  time  T  is 

TO  -  Y  l  r(t> dt  =  ¥  {  =  <B2» 

Y  l  'T-fiT  dt  =  4  [lQ(11'(t))]T.  " 

In  (1  -  F(T)) . 

A  system  consisting  of  two  independent  subsystems  in  parallel  where  the  subsystems 
have  exponential  durability  with  parameters  Xj  and  X2  is  an  example  of  a  system  which 
has  IFRA  durability  but  not  IFR  durability.  For  such  a  system, 


F(T)  = 

(1  -  e-X‘T) 

(1  -  e' 

-X2T)  5 

f(T)  = 

X,  e-XlT  (1 

-  e-X- 

r)  +  x2  e"XsT 

(1  -  e-X*T) , 

(B3) 

rlTl  — 

KI1 

_ 

e~XlT  -f  X2  e” 

X*T  -  (Xj  +  X2)  e“{X2  +  Xa)T 

rlA ;  — 

1  -  F(T) 

e_x*T  + 

e-X8T  _  e-(X»  +  X2)T 

T(T)  = 

~Y ln  [* 

-F(T) 

]=> 

(B4) 

31 


The  failure  rate  and  average  failure  rate  for  such  a  system  with  Xi  =  .3  and  X2  = 
.7  are  plotted  on  Figure  B.l. 


and 


so  that 


If  a  system  has  IFRA,  then  for  any  M  satisfying  the  condition  M  >  1,  T/M  <  T 


r(T/M)  =  <  -Minimi  =  r  (T) 


In  (1  -  F(T/M))m  >  ln(l  -  F(T)) 

°r 

[l  -  F(T/M)]  >  1  -  F(T) 

OF  M 

[r(T/M)]  >R(T). 


(B5) 


(B8) 


(B7) 


Equation  (B7)  can  be  used  to  construct  a  100(1  -  a)%  lower  confidence  limit  for 
R(T*)  under  the  IFRA  model,  where  T*  is  a  durability  life  requirement. 

If  N  units  are  tested  for  time  r  >  T*  and  r  units  survive  the  test,  then  we  can  con¬ 
struct  a  distribution  free  lower  confidence  limit  (frequentist  or  Bayesian)  for  R(T*) 
which  has  the  property 

P  |  R(r)  <  R(a,r)  J  =  a.  (B8) 


Since  t  >  T*  we  can  write  T*  =  T/M,  with  M  >  1.  It  then  follows  from  equa¬ 
tions  (B7)  and  (B8)  that 

P  |  R(t/M)m  <  R(a,r)  }  <  a  (B9) 

and 

P  {  R(r/M)  =  R(T*)  <  R(a,T*)1/M  J  <  a  ,  (BIO) 

that  is,  R(a,r)1/M  is  an  at  least  100(1  -  a)%  lower  confidence  limit  for  durability  relative 
to  durability  life  requirement  T*. 

Similar  arguments  can  be  developed  to  construct  confidence  limits  for  the  mean, 
median  and  other  quantiles  of  IFRA  systems.  Barlow  and  Proschan  (1981)  give  an 
excellent  treatment  of  this  subject. 
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Failure  Rate  arid  Failure  Rate  Average 


Figure  Bl.  Failure  rate  and  average  failure  rate  for  a  system  consisting 
of  two  independent  subsystems  with  exponential  durability 
in  parallel.  In  this  case,  both  subsystems  are  IFR  but 
the  system  is  IFRA  because  r(T)  is  increasing  even  though 
r(T)  is  decreasing  for  T  >  3.86. 
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APPENDIX  C.  THE  EXPONENTIAL  DISTRIBUTION 


The  exponential  distribution  is  one  of  the  most  useful  probability  distributions  for 
modeling  in  life  testing.  It  adequately  characterizes  many  realistic  processes  and.  its 
mathematical  simplicity  makes  it  easy  to  use.  It  also  has  the  property  that  it  bounds 
the  distributions  which  are  in  the  IFR  and  IFRA  families.  These  distributions  are  basic 
to  the  theory  of  durability  testing.  It  is  this  last  property  which  makes  it  important  to 
the  work  described  in  this  report. 

For  the  sake  of  completeness,  some  of  the  mathematical  properties  of  the  exponen¬ 
tial  distribution  are  developed  in  this  appendix.  Examples  which  demonstrate  the  use  of 
exponential  distribution  theory  are  also  given. 

The  user  of  the  procedures  developed  in  this  report  need  not  fret  with  the  details 
given  in  this  appendix  unless  he  is  motivated  by  scientific  curiosity.  The  remainder  of 
the  report  can  be  used  without  referring  to  this  appendix. 

A  random  variable,  T,  is  said  to  have  the  exponential  life  distribution  if  its  proba¬ 
bility  density  function  is  of  the  form 

f(T)#  =  4  e'T/* ;  T  >  0 , 0  >  0  .  (Cl) 

0 

It  follows  from  (Cl)  that  the  probability  that  a  unit  with  an  exponential  life  distri¬ 
bution  will  fail  at  or  before  time  T  is  given  by  the  distribution  function 

T  i 

F(T)  =  f  j-  e-T/'  dT  =  1  -  e-T/';  T  >0  .  (C2) 

o  0 

The  probability  that  the  life  of  such  a  unit  will  be  greater  than  T  is 

R(T)  =  1  -  F(T)  =  e-'1'/* .  (C3) 

In  this  report  we  often  refer  to  R(T)  as  the  durability  function. 

The  mean  of  the  exponential  distribution  is 

00  rp 

/  4  e"17'  dT  =  6  ;  (C4) 

0  " 

and  the  variance  is 

/  (T— )--  e-T/'  dT  =  &  .  (C5) 

o  0 

The  q’th  quantile,  Tq  ,  is  given  by  solving  for  Tq  in 


^  This  probability  density  function  is  sometimes  expressed  in  the  form  f(T)  —  X  e  We  chose  the  form  (Cl)  so  that  the 
expected  valne  of  T  would  be  6  rather  than  l/X. 


37 


p  {t<T,}  =q=l-e-T''\  (C6) 

which  has  solution  Tq  =  -0  In  (1  -  q)  =  -In  (1  -  q)#.  Thus,  the  median  of  the  distribu¬ 
tion  is 


T  5  =  -0  In  (.5)  =  (.6931)  0 ,  (C7) 

that  is,  the  median  is  less  than  the  mean  and  the  distribution  has  positive  skewness,  i.e., 
it  is  skewed  right  (  see  Figure  Cl). 

The  important  characterizations  of  durability  life  are  0,  the  average  durability  life; 
R(T),  the  probability  that  the  durability  life  of  a  randomly  selected  item  will  exceed 
time  T;  and  Tq,  the  q’th  quantile  of  the  durability  life  distribution.  Recall,  Tq  is  that 
life  which  will  be  exceeded  with  probability  1  -  q. 

Estimation  of  0,  R(T)  and  Tq  is  somewhat  of  a  problem  in  durability  testing.  In 
the  durability  life  problem  feasible  test  time  is  usually  short  relative  to  the  actual  dura¬ 
bility  life  of  an  item.  That  is,  if  a  vehicle  has  an  actual  durability  life  of  100,000  miles 
it  is  unlikely  that  it  would  be  practical  to  test  it  that  long.  Practicality  dictates  the 
selection  of  a  test  termination  (truncation)  time,  r.  N  items  are  tested  until  each  either 
operates  for  time  r  or  suffers  a  durability  failure.  If  C  items  fail,  the  total  test  time  for 
all  items  is 

A=  £  Ti  +  (N-C)r,  (C8) 

i=l 

where  Tj  is  the  failure  time  for  the  i’th  failed  item. 


The  likelihood  function  for  our  test  with  truncation  at  time  r  is 


Now 


so  that 


ln(L)  =  -  C  ln(0)  -  A/0 


d  i  n  \  -  ~C  A  _  o 
dO  D(  *  0  +  02  0 

and  the  maximum  likelihood  estimate  of  0  is 


(CO) 


0  =  A/C.  (CIO) 

where  A  is  given  in  relationship  (C8).  The  above  estimate  of  0  can  be  used  to  estimate 
other  characteristics  of  durability  life  (see  Table  C-l). 
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EXPONENTIAL  DISTRIBUTION 


A 

We  have  found  that  0  ,  as  defined  above,  is  approximately  distributed  as  a  chi- 
square  variate  with  2C  +  1  degrees  of  freedom.*  This  approximation  is  used  to  give  the 
lower  confidence  limits  in  Table  C-l. 

Bartholomew  (1963)  derived  the  exact  distribution  of  0  when  C^O.  A  computer 
program  which  uses  his  results  to  construct  exact  lower  confidence  limits  for  9  is  given 
in  Appendix  E. 

Clearly,  the  maximum  likelihood  estimate  of  0  is  not  defined  when  C  =  0  (see 
Equation  CIO). 

TABLE  C-l.  Estimates  of  durability  distribution  properties  for  time 
truncated  test  with  one  or  more  failures. 


PROPERTY 

Mean  (Ave) 

Durability 

(POINT)  ESTIMATE# 

•=iA 

APPROX.  L.C.L.  |100(l-a)%] 

e  _  2  A 

L  X2i-a  (2C  4-  1) 

q’th  quantile 

of  durability  life 

H 

II 

SL 

1 

q,  _  -2  A  ln(l  -  q) 

U  “  X2l-a(2C  +  1) 

Median 

durability  life 

T  5  =  6  ln(2) 

^  _  2  A  ln(2) 

X2i-a  (2C  4-  1) 

Durability  for  test 

time  t 

R(r)  =  e-r/* 

D,  x  f-(r)X2i-a  (2C  +  1)  ] 

R(a,r)  =  exp  - — - 

Durability  for 

mission  time 

T*  =  r/M 

R(r/M)  =  e-f/M*  =  [r(T*)]  */M 

R(a,T*)  =  |^R(a,T*)] 

#  Test  truncated  at  time  r  with  one  or  more  durability  failures.  Durability  mission  time 
is  T*  =  r/M,  C  =  number  of  durability  failures  and  A  =  total  test  time  = 

E  Tj  4-  (N-C)  r. 
isi _ 

There  are  many  approximations  for  this  confidence  limit.  The  one  which  is  most  common  involves  the  use  of  a  chi-square 
variate  with  2C+2  degrees  of  freedom.  Preliminary  results  from  our  study  of  the  problem  suggest  that  a  chi-square  with  2C+1  d.f., 
which  has  been  considered  by  others,  may  be  the  best  chi-square  approximation.  At  this  time  we  recommend  using  the  chi-square 
approximation  with  2C+1  d.f.  because  it  gave  confidence  limits  which  are  closest  to  the  exact  limits  in  a  recent  monte-carlo  study. 
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However,  we  can  construct  lower  confidence  limits  by  combining  distribution  free  pro¬ 
cedures  (see  Appendix  A)  with  properties  of  the  exponential  distribution. 


First,  we  note  that  a  100(1  -  a)%  distribution  free  lower  confidence  limit  for  R(r)  is 
given  by  the  solution  of  the  equation 

E  (?)  F R(«.  r)  1 '  (l  -R(«,  r)  1  =  a.  (C1I) 

i=N-C  L  J  L  J 

Thus,  when  C  =  0, 

R(a,  r)  =  a!/N  .  (C12) 

We  next  take  another  approach  to  this  problem.  Assume  that  0,  if  it  existed,  has 
the  approximate  chi-square  distribution  with  2C,  rather  than  2C+1  degrees  of  freedom 
and  that,  if  we  had  continued  testing  for  time  r  >  +5,  5  =  0,  there  would  have  been 
one  durability  failure  at  that  time.  It  is  easy  to  show  that  x2i_<»  (2)  =  -2  In  a.  In  this 
case  A  =  N  r.  From  Table  C-l,  we  see  that 

R(a,r)  =  e-+x*U*)/u  (C13) 

_  £-f{-2Ina)/2Nr 

_  elnar1/N 

=  o‘/N  , 

which  is  the  same  as  the  distribution  free  lower  confidence  limit  for  R(t). 

Using  this  approach,  i.e.,  substituting  -2  In  a  for  X2i-a  (2C  +  1)  in  Table  C-l,  leads 
to  the  lower  confidence  limits  for  the  case  C  =  0  which  are  given  in  Table  C-2. 
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TABLE  C-2.  Lower  confidence  limits  for  durability  distribution  properties 
for  a  time  truncated  test  with  zero  failures. 


PROPERTY 

APPROXIMATE  100(l-o)%  L.C.L. 

Mean  (Ave.) 
Durability 

0L  =  A/lno 

q’th  quantile 
of  reliability  life 

TqL  =  A  ln  i1  ~  q)/10  a 

Median  durability 
life 

T. 5,  =  -  A  la(2)/ln  a 

Durability  for  test 
time  t 

R(a,r)  =  a1/N 

Durability  for 
mission  time 

T*  =  t/M 

1  ' 
R(a,T*)  =  a™ 
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APPENDIX  D.  BILLINGS’  MODERATELY  DISTRIBUTION  FREE  PROCEDURE 


Billings  (1967)  placed  the  requirement 

Tf(T)  >  F(T) ;  0  <  T  <  T0  (Dl) 

on  the  durability  probability  distribution  and  density  functions.  A  pictorial  representa¬ 
tion  of  this  assumption  is  given  in  Figure  Dl. 


If  (Dl)  is  satisfied,  then  it  follows  that,  on  the  interval  (0,  T0] 


f(T) 


> 


and  for  any  number  M  >  1, 


so  that 


or,  equivalently, 


which  leads  to 


F(T)  -  T 


7  >  7  4^ 


T  F(T)  t  X 
In  F(MT)  -  In  F(T)  >  In  (MT)  -  ln(T) 

In 


F(MT) 

>  In 

MT 

F(T) 

T 

C 

F(MT)  >  F(T) . 


(D2) 


(1)3 


(Dl) 


(D5) 


(DO) 


Billings  used  equation  (D6)  to  construct  confidence  bounds  on  durability  life  for  a 
stated  durability  life  requirement,  T*. 


Suppose  each  of  N  units  is  tested  until  it  fails  or  for  time  MT*,  whichever  comes 
first.  The  probability  that  any  particular  unit  will  experience  a  durability  failure  during 
the  test  is  F(MT*).  If  Y  units  fail  during  the  test  then  a  100(1  -  oJ'Y'  upper  confidence 
limit  for  F(MT*)  is  given  by  the  solution  of  the  equation 

£  C|  F„(MT*)X  [l-FL  (MT.)l  =  a  (D7) 

X— 0  "  L  J 

for  Fv  (MT*). 

If  Fu  (MT*)  is  the  solution  of  (D7)  then  we  can  make  the  probability  statement 


Since  ^4F(T*)  <  F(MT*)  it  follows  that 

P  (MF(T*)  >  F,j  (MT*)  |  <  a  (DO) 
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Figure  D1  .  Illustration  of  a  durability  distribution 

for  which  Tf(T)  >_  F(T)  for  all  T  satisfying 
0  <  T  Tq.  Here  T  f(T^)  is  the  area  under 
the  rectangle  bounded  by  the  axes  and  the 
dashed  lines  and  F(T^)  is  the  area  in  the 
hatched  region  under  the  curve. 


so  that 


=  P  Jf(T*)  >  Fu(MT*)/m} 

=  P  1 1  -  F(T*)  <  1  -  Fu  (MT*)/M  J 

P  Jr(T*)  <  1  -  Fu  (MT*)/M  J  <  a  .  (DIO) 

Thus,  1  -  Fu  (MT*)/M  =  R(aT*)  is  at  least  a  100(1  -  a)%  lower  confidence  limit 
for  R(T*),  which  is  the  probability  that  the  durability  life  is  a  randomly  selected  unit 
that  will  exceed  the  requirement,  T*. 

For  the  special  case  in  which  no  units  fail  during  the  test  time  MT*,  Y  =  0  and 
equation  (D7)  had  the  form 

£l-Fu(MT*)j  =a  (DU) 

so  that 

Fu  (MT*)  =  1  -  a‘/N  .  (D12) 

It  follows  from  equation  (DIO)  that  1  -  (1  -  ax/N)/M  is  an  at  least  100(1  -  a)% 
lower  confidence  limit  for  R(T*). 

Billings  (1968)  uses  equation  (D7)  to  refine  the  probability  statement  made  in  equa¬ 
tion  (DIO)  when  Y  >  0.  This  work  leads  to  the  relationship 

P  (r(T.)  <  1  -  -i-  F„  (MT.) }  <  a  ,  (D13) 

where  k  >  1  depends  on  N  and  a.  The  details  of  this  work,  which  are  mathematically 
tedious,  are  not  given  here. 
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APPENDIX  E.  LISTINGS  OF  COMPUTER  PROGRAMS 


49 


ooooooooop 


APPENDIX  E.  LISTINGS  OF  COMPUTER  PROGRAMS 


THIS  PROGRAM  FINDS  THE  LOWER  CONFIDENCE  LIMIT  (THETA) 
FOR  A  DURABILITY  STUDY  GIVEN  : 

ALPHA 

SAMPLE  SIZE  (N) 

TERMINATION  TIME  (T) 

AND  THETA(NOT) 

ALPHA  =  PROB  (  THETA  .GE.  THETA(NOT) ) 


PROGRAM  CONF  ( INPUT, OUTPUT,TAPE5=INPUT,TAPE6=OUTPUT,TAPE8) 
DIMENSION  IRAY(6),THETA(1) 

REAL  EPS,EPS2,ETA 
INTEGER  NSIG,NO,ITMAX,IER 
COMMON  ALPHA, N,T,THENOT 
EXTERNAL  F 
DATA  IRAY/6*  -0/ 

C - SET  PRINT  LIMIT  TO  ZERO 

IRAY(4)  =  0 

t - ALPHA  =  PR(THETA  HAT  .GE.  THETA  NOT)  = 

C - 1/  (1-EXP(-N*T/THETA))  [SUM  FROM  K=1  TO  N[(COMB(N,K)  • 

C - [SUM  FROM  1=0  TO  K[(COMB(K,I)  *  (-1)**I  * 

C - EXP(-T/THETA(N-K+lj)  * 

C - INTEGRAL  FROM  X  TO  INF  OF  P(CHI**2  W  2K  DF)  D(CHI**2)  2K 

WRITE  (6,100) 

100  FORMAT  (’ ALPHA  LEVEL’) 

READ  (5,*)  ALPHA 
WRITE  (6,110) 

110  FORMAT  (’ SAMPLE  SIZE’) 

READ  (5,*)  N 
WRITE  (6,120) 

120  FORMAT  (’ TRUNCATION  TIME’) 

READ  (5,*)  T 
WRITE  (6,130) 

130  FORMAT  (’ THETA(NOT)’) 

READ  (5,*)  THENOT 


EPS  =0.0 
NSIG  =  4 

XL  =  THENOT/20 
XR  =10 
ITMAX=  100 

C - CALL  SYSTEMC  TO  INHIBIT  PRINTING  OF  ERROR  115 

C - AND  START  ERROR  SUMMARY  ACCUMULATION 
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CALL  SYSTEMC(115,IRAY) 

CALL  ZFALSE  (F,EPS,NSIG,XL,XR,XAPP,ITMAX,IER) 

WRITE  (6,600)  XAPP, ALPHA, N,T,THENOT 
600  FORMAT  (’  THE  LOWER  CONFIDENCE  LIMIT  —  ’.F12.5// 
•’  FOR  ALPHA  =  ’,F6.4/ 

•’  N  -  ’,16/ 

•’  T  *=  \F12.6/ 

*’  THETA(NOT)  =  ’.F12.6) 


STOP 

END 


C - 

REAL  FUNCTION  F(THETA) 

REAL  THETA 

REAL  A, B,T,X,DF,P, IERR,FACT1,FACT2,SUM1,SUM2, TEETHE, VALUE 
REAL  PRODl,PROD2,Z, ALPHA, THENOT, TOTAL, EPFJEPF2, ETA 
INTEGER  N,K,I,NN,KK,Il,NSIG,ITMAX,NO 
COMMON  ALPHA, N,T,THENOT 

A  =  1  /  (1  -  EXP(-N*T/THETA)  ) 

TEETHE  =  -T  /  THETA 

C - SUM2  =  [SUM  FROM  K=1  TO  N)(COMB(N,K)  * 

C - [SUM  FROM  1=0  TO  K](COMB(K,I)  *  (-1)**I  * 

C - EXP(-T/THETA(N-K+I)) 

SUM2  =  0 
FACT2  =  1.0 
PROD2  —  0 
NN  =  N 
DO  10  K  =  1,N 
FACT2  =  FACT2  *  NN/K 

SUM1  =  0 
PROD1  =  0 
KK  =  K 
DO  20  I  =  0,K 
IF  (I  .EQ.  0)  THEN 
FACT1 = 1 
ELSE 

FACT1  =  FACT1  *  KK/I 
KK  =  KK -  1 
END  IF 

PROD1  =  FACT1  *  (-1)**I  *  EXP  (TEETHE  *  (N-K+I)) 

C— Z  =  l-(INT(FROM  X  TO  INF)  OF  P(CHI**2  W  DF=2K)  D(CHI**2)  DK) 


DF  =  2*K 

VALUE  =  THENOT  -  T/K  *  (N-K+I) 
IF  (VALUE  .LE.  0)  THEN 
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Z-  1. 

ELSE 

IF  (THETA  .LT.  0)  THEN 
Z  =  1 
GO  TO  50 
ELSE 
END  IF 

X  =  2*K/THETA  *  (VALUE) 
CALL  MDCH(X,DF,P,IERR) 

Z  —  1-P 

C  WRITE(6,*)  Z 

END  IF 

50  PROD1  «  PROD1  *  Z 
SUM1  =  SUM1  +  PROD1 
20  CONTINUE 

PROD2  =  FACT2  *  SUMl 
SUM2  =  SUM2  +  PROD2 
NN  —  NN-1 


10  CONTINUE 

F  =  (A  *  SUM2)  -  ALPHA 

RETURN 

END 


PROGRAM  USED  TO  GENERATE  TABLE  F-l 


C  THIS  PROGRAM  GIVES  A  TABLE  OF  THE  SAMPLING  PLAN  FOR  THE 
C  DURABILITY  OF  IFRA  DISTRIBUTIONS 
C 

C  INPUTS  ARE: 

C  DESIRED  CONFIDENCE  LEVEL  (DCL)  AND  NUMBER  OF  FAILURES  (C) 

PROGRAM  BARLOW  (INPUT, OUTPUT, TAPE5=  INPUT, TAPE6=OUTPUT,TAPE8  TAPE9) 
DIMENSION  NTOTAL(20,20),TRUECL(20,20),TIME(16),DUR(6) 

INTEGER  C 

REAL  P,DCL,PDF,CDF,TCDF,R,Tl 
C  N TOTAL  =  TOTAL  SAMPLE  NUMBER 

C  TRUECL  =  TRUE  CONFIDENCE  LEVEL 
C  TIME  =  THE  MULTIPLE  OF  TIME  THAT  IS  MULTIPLIED  BY 
C  DURABILITY  CRITERION  TO  GET  TOTAL  TEST  TIME 

C  DUR  =  DURABILITY  CRITERION 

C  INPUT  VALUES  FOR  TIME  AND  DURABILITY  CRITERION: 

DATA  TIME/l. 0,1. 1,1. 2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 

+  2.0, 2.2, 2.4, 2.6, 2.8, 3.0/ 

DATA  DUR/. 50, .60, .70,. 80,. 90,. 95/ 

WRITE  (6,100) 

100  FORMAT  (’DESIRED  CONFIDENCE  LEVEL  ’) 

READ  (5,*)  DCL 
WRITE  (6,110) 

110  FORMAT  (  NO.  OF  FAILURES  ’) 

READ  (5,*)  C 

C  TABLE  SET-UP 

WRITE  (9,120)  DCL,C 

120  FORMAT  (T,25X, ’SAMPLING  PLAN  FOR  THE  DURABILITY  OF  IFRA' 

+’  DISTRIBUTIONS’, //,37X,’  DESIRED  CONFIDENCE  LEVEL  =  ’, 

+F3.2,/,37X,’  NUMBER  OF  FAILURES  (C)  =  ’,11,//) 

WRITE  (9,130) 

130  FORMAT  (1X,98(’-’)) 

WRITE  (9,140) 

140  FORMAT  (’  :’,10X,’:’,30X, ’DURABILITY  ’,44X,Y) 

WRITE  (9,203) 

WRITE  (9,150)  (DUR(M),M=1,6) 

150  FORMAT  (’  :’,10X,’:’,5X,5(F3.2,11X),F3.2,7X,’:’) 

WRITE  (9,190) 

190  FORMAT  (’:  M*  :’,85(’-’),’:’) 

WRITE  (9,170) 

WRITE  (9,180) 

170  FORMAT  (’:  :  ’,5(’  N  TRUE  ’),’  N  TRUE  :’) 

180  FORMAT  (’  :  :  ’,5(’  C.L.  ’),’  C.L.  :’) 

WRITE  (9,130) 

WRITE  (9,203) 


54 


C  FOR  EACH  MULTIPLE  OF  TIME  AND  EACH  DURABILITY  CRITERION, 

C  THE  NUMBER  OF  TEST  ITEMS  (SAMPLE  SIZE)  AND  TRUE  CONFIDENCE 
C  LEVEL  ARE  CALCULATED. 

C 

PDF  =  0 
CDF  =  0 
PPP  =  1  -  DCL 

DO  10  I  =  1,16 
T1  =  TIME(I) 

N  =  C  +  1 

DO  20  J  =  1,6 
R  =  DUR(J) 

P  =  1  -  (R**T1) 

55  CALL  BINOM  (PDF,CDF,N,P,C) 

IF  (CDF  .LE.  PPP)  THEN 
TCDF  =  1-CDF 
NTOTAL(I,J)  =  N 
TRUECL(I,J)  =  TCDF 
ELSE 
N  —  N+l 
GO  TO  55 
END  IF 

20  CONTINUE 
10  CONTINUE 

C  THE  NUMBER  OF  UNITS  AND  THE  TRUE  CONFIDENCE  LEVEL 
C  ARE  PRINTED  IN  THE  TABLE 

I  =  1-1 
J  =  J-l 
DO  40  II  =  1,1 

WRITE  (9,200)  TIME(II),(NTOTAL(II,JJ),TRUECL(II,JJ),JJ=l,J) 

200  FORMAT  (’  :’,2X,F4.1,4X,’:’,1X,6(I4,2X,F4.2,4X),’:’) 

WRITE  (9,203) 

203  FORMAT  (’  :’,10X,’:’,85X,’:’) 

40  CONTINUE 
WRITE  (9,203) 

WRITE  (9,205) 

205  FORMAT  (IX, 98(’-’)) 

WRITE  (9,210) 

210  FORMAT  (IX,/,’  *  M  =  TEST  TIME/DURABILITY  CRITERION’/////) 
STOP 
END 


C  THIS  SUBROUTINE  FINDS  THE  PROBABILITY  DENSITY  FUNCTION 
C  AND  CUMMULATIVE  DENSITY  FUNCTION  GIVEN  THE  SAMPLE  SIZE 
C  (N),  THE  P-VALUE,  AND  THE  NUMBER  OF  FAILURES  (C). 

SUBROUTINE  BINOM  (PDF,CDF,N,P,C) 

REAL  PDF, CDF, P 
INTEGER  C 
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PDF  =  (1-P)**N 
CDF  =  PDF 
IF  (C  .NE.  0)  THEN 
PI  =  P/(l-P) 

DO  30  K  =  1,C 

PDF  =  PDF  *  (N-K+l)  •  Pl/K 
CDF  =  CDF  +  PDF 
30  CONTINUE 
ELSE 
END  IF 


RETURN 

END 


PROGRAM  USED  TO  GENERATE  TABLE  G-l 


C  THIS  PROGRAM  CREATES  A  FILE  THAT  WILL  LIST  A  TABLE  OF 
C  A  SAMPLING  PLAN  FOR  THE  MEAN  OF  IFRA  DISTRIBUTIONS. 

C  NTOTAL  =  NUMBER  OF  TEST  ITEMS  (  SAMPLE  NUMBER  ) 

C  TRUECL  =  TRUE  CONFIDENCE  LEVEL 

C  TIME  =  MULTIPLE  OF  TIME  THAT  THE  DURABILITY  CRITERION 
C  IS  MULTIPLIED  BY  TO  OBTAIN  THE  TOTAL  TEST  TIME 

C  FAILURE  =  C  =  NUMBER  OF  FAILURES  ALLOWED  IN  A  SAMPLE 

C  C  RANGES  FROM  0  TO  5 

PROGRAM  BARLOW  (INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE8,TAPE9) 
DIMENSION  NTOT AL(20, 20), TRUECL(20, 20), TIME(15), FAILURES  1) 

INTEGER  C, FAILURE 

REAL  P,DCL,PDF,CDF,TCDF,R,T1,A,L 

DATA  TIME/1. 1,1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 

+  2.0, 2.2, 2.4, 2.6, 2.8, 3.0/ 

DATA  FAILURE/0, 1,2, 3, 4, 5, 6, 7, 8, 9, 10/ 

C  INPUT  THE  DESIRED  CONFIDENCE  LEVEL 

WRITE  (6,100) 

100  FORMAT  (’DESIRED  CONFIDENCE  LEVEL  ’) 

READ  (5,*)  DCL 

C  WRITE  THE  TITLE  AND  COLUMN  HEADINGS 
WRITE  (9,120)  DCL 

120  FORMAT  (T,25X, ’SAMPLING  PLAN  FOR  THE  MEAN  OF  IFRA’, 

+’  DISTRIBUTIONS’, //,37X,’  DESIRED  CONFIDENCE  LEVEL  =  ’, 

+F3.2 ,//) 

WRITE  (9,130) 

130  FORMAT  (1X,98(’-’)) 

WRITE  (9,140)  (FAILURE(M),M=1,6) 

140  FORMAT  (’  :’,10X,’:’,2X,’C  —  ’,I2,12X,4(I2,12X),I2,7X,’:’) 

WRITE  (9,150) 

150  FORMAT  (’  :  F(TIME)*  :’,85(’-’),’:’) 

WRITE  (9,170) 

WRITE  (9,180) 

170  FORMAT  (’  :  :  ’,5(’TEST  TRUE  ’),’TEST  TRUE  :’) 

180  FORMAT  (’  :  :  ’,6(’ITEMS  C.L.  ’),’:’) 

WRITE  (9,130) 

WRITE  (9,190) 

190  FORMAT  (’  :’,10X,’:’,85X,’:’) 


C  PDF  =  PROBABILITY  DENSITY  FUNCTION 
C  CDF  =  CUMMULATIVE  DENSITY  FUNCTION 


PDF  =  0 
CDF  =  0 
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DO  10  I  =  1,15 
T1  =  TIME(I) 

A  =  1  -  DCL 
DO  20  L  =  A,l,.l 
F  =  1-L-EXP(-L*T1) 

IF  (F  .GT.  0)  THEN 
GO  TO  35 
ELSE 
A  =  L  -  .1 
B  =  L 
GO  TO  25 
END  IF 

35  CONTINUE 
20  CONTINUE 

C  FIND  W 

25  CALL  W1(A,B,T1,W) 

C  FOR  EACH  VALUE  OF  THE  NUMBER  OF  FAILURES,  C,  FIND 
C  THE  TRUE  CONFIDENCE  LEVEL,  TCDF,  SUCH  THAT  TCDF 
C  IS  GREATER  THAN  OR  EQUAL  TO  THE  DESIRED  CONFIDENCE 
C  LEVEL,  DCL. 

DO  30  J  =  1,6 
C  =  FAILURE(J) 

P  =  1-EXP(-W  *  Tl) 

IF  (C  .EQ.  0)  N  =  1 

360  CALL  BINOM  (PDF,CDF,N,P,C) 

IF  (CDF  .LE.  (1-DCL))  THEN 
TCDF  =  1  -  CDF 
NTOTAL(I.J)  =  N 
TRUECL(I,J)  =  TCDF 
ELSE 

N  =  N  +  1 
GO  TO  360 
END  IF 

30  CONTINUE 
10  CONTINUE 

C  WRITE  THE  TABLE  INPUT:  TOTAL  TIME/DURABILITY  REQUIREMENT, 
C  TEST  ITEMS  (SAMPLE  SIZE),  AND  TRUE  CONFIDENCE  LEVEL. 

I  =  1-1 
J=  J-l 
DO  40  II  =  1,1 

WRITE  (9,200)  TIME(II),(NTOTAL(II,JJ),TRUECL(II,JJ),JJ=l,J) 

200  FORMAT  (’  I’^XJ^.l^X/i’.lX.efM^XJ^.^X),’:’) 

WRITE  (9,190) 

40  CONTINUE 
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WRITE(9,190) 

WRITE  (9,210) 

210  FORMAT  (1X,98(’-’)) 

WRITE  (9,220) 

220  FORMAT  (IX,/,’  *  F(TIME)  =  TEST  TIME/DURABILITY  CRITERION’//) 
STOP 
END 


C  COMPUTES  W.  INPUT  VALUES  ARE  A,B,T1. 

SUBROUTINE  W1  (A,B,T1,W) 

REAL  A,B,T1,W 

IF  (ABS(1-A-EXP(-A*T1))  .GT.  IE-09  )  THEN 
IF  (  ABS(1-B-EXP(-B*T1))  .GT.  IE-09  )  THEN 
350  D  =  (A+B)/2 

IF  (ABS(A-B)  .LT.  1E-09)  THEN 
W  =  A 
ELSE 

F  =  1-D  -  EXP(-D*T1) 

IF  (  ABS(F)  .LT.  IE-09)  THEN 
W  =  D 
ELSE 

IF  (  F  .GT.  0  )  THEN 
A  =  D 
ELSE 
B  =  D 
END  IF 
GO  TO  350 
END  IF 
END  IF 
ELSE 
W  =  B 
END  IF 
ELSE 
W  =  A 
END  IF 
RETURN 
END 

C  COMPUTES  THE  CUMULATIVE  DISTRIBUTION  FUNCTION  AND 
C  PROBABILITY  DENSITY  FUNCTION  USING  THE  BINOMIAL. 

C  INPUT  VALUES  ARE  N,  P,  C. 

SUBROUTINE  BINOM  (PDF,CDF,N,P,C) 

REAL  PDF, CDF, P 
INTEGER  C 
PDF  =  (1-P)**N 
CDF  —  PDF 
IF  (C  .NE.  0)  THEN 
PI  =  P/(l-P) 

DO  50  K  =  1,C 

PDF  =  PDF  *  (N-K+l)  *  Pl/K 
CDF  —  CDF  +  PDF 
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50  CONTINUE 
ELSE 
END  IF 
RETURN 
END 

THIS  IS  PAGE  TWO  OF  THE  TABLE  FOR  ’SAMPLING  PLAN  FOR 
THE  MEAN  OF  IFRA  DISTRIBUTIONS’ 

C  RANGES  FROM  6  TO  10 


PROGRAM  BARLOW  (INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE8,TAPE9) 
DIMENSION  NTOTAL(20,20),TRUECL(20,20),TIME(15),FAILURE(11) 

INTEGER  C, FAILURE 

REAL  P,DCL,PDF,CDF,TCDF,R,T1,A,L 

DATA  TIME/l.  1,1. 2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 

+  2.0, 2.2, 2.4, 2.6, 2.8, 3.0/ 

DATA  FAILURE/0,1,2,3,4,5,6,7,8,9,10/ 

WRITE  (6,100) 

0  FORMAT  (’DESIRED  CONFIDENCE  LEVEL  ’) 

READ  (5,*)  DCL 


WRITE  (9,120)  DCL 

120  FORMAT  (’l’,22X, ’SAMPLING  PLAN  FOR  THE  MEAN  OF  IFRA’, 
+’  DISTRIBUTIONS’, //,31X,’  DESIRED  CONFIDENCE  LEVEL  =  ’, 
+F3.2,//) 

WRITE  (9,130) 

130  FORMAT  (1X,84(’-’)) 

WRITE  (9,140)  (FAILURE(M),M=7,11) 

140  FORMAT  (’  :’,10X,’:’,2X,’C  =  ’,I2,12X,3(I2,12X),I2,7X,’:’) 

WRITE  (9,150) 

150  FORMAT  (’  :  F(TIME)*  :’,71(’-’),’:’) 

WRITE  (9,170) 

WRITE  (9,180) 

170  FORMAT  (’:  :  ’,4(’TEST  TRUE  ’),’TEST  TRUE  :’) 

180  FORMAT  (’  :  :  ’,5(’ITEMS  C.L.  ’),’:’) 

WRITE  (9,130) 

WRITE  (9,190) 

190  FORMAT  (’  :’,10X,’:’,71X,’:’) 


PDF  =  0 
CDF  =  0 


DO  10  I  =  1,15 
T1  =  TIME(I) 

A  = 1  - DCL 
DO  20  L  =  A,l,.l 
F  =  1-L-EXP(-L*T1) 
IF  (F  .GT.  0)  THEN 
GO  TO  20 
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ELSE 
A  =  L  -  .1 
B  —  L 
GO  TO  25 
END  IF 

20  CONTINUE 

25  CALL  W1(A,B,T1,W) 

DO  30  J  =  7,11 
C  =  FAILURE(J) 

P  =  1-EXP  (-W  *  Tl) 

IF  (C  .EQ.  6)  N  =  1 

360  CALL  BINOM  (PDF,CDF,N,P,C) 

IF  (CDF  .LE.  (1-DCL))  THEN 
TCDF  =  1  -  CDF 
NTOTAL(I,J)  =  N 
TRUECL(I,J)  =  TCDF 
ELSE 

N  =  N  +  1 
GO  TO  360 
END  IF 

30  CONTINUE 
10  CONTINUE 


I  =  1-1 
J  =  J-l 
DO  40  II  =  1,1 

WRITE  (9,200)  TIME(II),(NTOTAL(tI,JJ),TRUECL(II,JJ),JJ=7,J) 

200  FORMAT  (’  :’,2X,F4.1,4X,’:’,1X,5(I4,2X,F4.2,4X),V) 

WRITE  (9,190) 

40  CONTINUE 

WRITE  (9,190) 

WRITE  (9,210) 

210  FORMAT  (1X,84(’-’)) 

WRITE  (9,220) 

220  FORMAT  (IX,/,’  *  F(TIME)  =  TEST  TIME/DURABILITY  CRITERION’//) 
STOP 
END 

SUBROUTINE  W1(A,B,T1,W) 

REAL  A,B,T1,W 

IF  (ABS(1-A-EXP(-A*T1))  .GT.  IE-09  )  THEN 
IF  (  ABS(1-B-EXP(-B*T  1 ))  .GT.  IE-09  )  THEN 
350  D  =  (A+B)/2 

IF  (ABS(A-B)  .LT.  IE-09)  THEN 
W  =  A 
ELSE 

F  =  1-D  -  EXP(-D*T1) 

IF  (  ABS(F)  .LT.  IE-09)  THEN 
W  =  D 
ELSE 
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IF  (  F  .GT.  0  )  THEN 
A  =»  D 
ELSE 
B  =  D 
END  IF 
GO  TO  350 
END  IF 
END  IF 
ELSE 
W  =  B 
END  IF 
ELSE 
W  =  A 
END  IF 
RETURN 
END 

SUBROUTINE  BINOM  (PDF,CDF,N,P,C) 
REAL  PDF, CDF, P 
INTEGER  C 
PDF  —  (1-P)**N 
CDF  =  PDF 
IF  (C  .NE.  0)  THEN 
PI  =  P/(l-P) 

DO  50  K  =  I,C 

PDF  =  PDF  *  (N-K+l)  *  Pl/K 
CDF  —  CDF  +  PDF 
50  CONTINUE 
ELSE 
END  IF 
RETURN 
END 
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PROGRAM  USED  TO  GENERATE  TABLE  H-l 


PROGRAM  BARLOW  (INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE8,TAPE9) 
DIMENSION  TIME(20,20),  CONF(9),DURAB(8),CC(2),NN(9) 

REAL  P, CDF, M,R 
INTEGER  N,C,NN,CC 

C  THIS  PROGRAM  GIVES  A  TABLE  OF  A  SAMPLING  PLAN  FOR 
C  THE  MISSION  TIME  FOR  IFRA  DISTRIBUTIONS. 

C  THIS  IS  DONE  SIMPLY  BY  INPUTTING  SAMPLE  SIZE  (N)  AND 
C  NUMBER  OF  FAILURES  (C). 

C  TIME  =  A  MULTIPLE  OF  TIME  THAT  WILL  BE  CALCULATED  SO 
C  TIME  MULTIPLIED  BY  DURABILITY  CRITERION  EQUALS 

C  TOTAL  TEST  TIME 

C  CONF  =  REQUIRED  CONFIDENCE  LEVEL 
C  DURAB  =  R  =  DURABILITY  CRITERION 
C  C  =  NUMBER  OF  FAILURES 

C  N  =  SAMPLE  SIZE 

DATA  DURAB/. 50, .60, .70, .75, .80, .85, .90, .95/ 

DATA  CONF /.50,. 60,. 70,. 75, .80, .85,. 90,. 95, .99/ 

DATA  CC/0,1/ 

DATA  NN/2,3,4,5,6,7,8,9,10/ 


C  TITLE  AND  TABLE  FORMAT.  DURABILITY  AND  CONFIDENCE  LEVELS 
C  ARE  ON  THE  X  AND  Y  AXES  RESPECTIVELY. 

DO  70  KK  =  1,2 
C  =  CC(KK) 

DO  80  LL  =  1,9 
N  =  NN(LL) 

WRITE  (9,120)  N,C 

120  FORMAT  (T, ’TABLE  II.  (CONT”D)’,10X, ’SAMPLING  PLAN  FOR’, 

+’  M  FOR  IFRA’, 

+’  DISTRIBUTIONS’, //,38X,’  NUMBER  OF  UNITS  (N)  =  ’, 

+I2,/,38X,’  NUMBER  OF  FAILURES  (C)  —  ’,11,//) 

WRITE  (9,130) 

130  FORMAT  (1X,98(’-’)) 

WRITE  (9,140) 

140  FORMAT  (’  :’,10X,’:’,31X, ’DURABILITY’, 44X,’:’) 

WRITE  (9,150) 

150  FORMAT  (’  :’,10X,’:’,85X,’:’) 

WRITE  (9,160)  (DURAB(L),L=  1,8) 

160  FORMAT  (’  CONF.  5X,7(F3.2,7X),F3.2,7X,’:’) 

WRITE  (9,150) 

WRITE  (9,130) 

WRITE  (9,150) 


FOR  EACH  CONFIDENCE  AND  DURABILITY  REQUIREMENT,  TIME  IS 
CALCULATED  (THAT  IS  THE  MULTIPLE  OF  DURABILITY  CRITERION 
WHOSE  PRODUCT  IS  TOTAL  TIME). 


63 


DO  10  I  —  1,9 
CDF  —  CONF(I) 
K  =  0 


DO  20  J  =  1,8 
R  =  DURAB(J) 

K  =  K  +  1 

CALL  BINOM  (CDF,N,R,C,M,K) 
IF  (  M  .LT.  1.0  )  THEN 
M  =  M  *  10**8 
ELSE 
END  IF 

TIME  (1,3)  =  M 
20  CONTINUE 
10  CONTINUE 

C  TABLE  VALUES  (TIME)  PRINTED. 


I  =  1-1 
J  =  J-l 
DO  40  II  =  1,1 

WRITE  (9,200)  CONF(II),(TIME(II,JJ),JJ=l,J) 

200  FORMAT  (’  :\3X,F3.2,4X,’:’,3X,7(F6.3,4X),F6.3,6X,’:’) 

WRITE  (9,150) 

40  CONTINUE 
WRITE  (9,150) 

WRITE  (9,210) 

210  FORMAT  (IX, 98(’-’)) 

WRITE  (9,220) 

220  FORMAT  (//’  ******  =  M  <  1’, 

+/’  RECOMMENDATION:  1  <  M  <  3  ’,/////) 

80  CONTINUE 
70  CONTINUE 
STOP 
END 


C 

C  CALCULATES  M  —  TIME  GIVEN  CUMMULATIVE  DIST’N  FUNCTION, 
C  SAMPLE  SIZE  (N),  DURABILITY  (R),  AND  #  OF  FAILURES  (C). 

SUBROUTINE  BINOM  (CDF,N,R,C,M,K) 

REAL  CDF,R,M,P,F,I,ALPHA,RO 
INTEGER  C,N 
ALPHA  =  1  -  CDF 
XN  =  N 

IF  (C  .EQ.  0)  THEN 

M  =  (ALOG(ALPHA))/  (  XN  *  ALOG(R)) 

ELSE 

IF  (C  .EQ.  1)  THEN 
IF  (K  .EQ.  1)  THEN 
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50 


RO  =  0 

I  —  .1 

RO  —  RO  +  I 

TEST  =  RO**XN  +  XN  *  (1-  RO)  *  RO**(XN-l) 
IF  (  ABS(TEST  -  ALPHA)  .LT.  .000001  )  THEN 
GO  TO  60 
ELSE 

IF  (  TEST  .LT.  ALPHA)  THEN 
RO  =  RO  +  I 
ELSE 

RO  =  RO  -  I 
I  =  I  /  10. 

RO  —  RO  +  I 
END  IF 
END  IF 
GO  TO  50 
ELSE 
END  IF 
ELSE 
END  IF 
60  CONTINUE 

M  —  (ALOG(RO))  /  (ALOG(R)) 

END  IF 

RETURN 

END 
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APPENDIX  F.  TABLE  F-l.  SAMPLING  PLAN  FOR  THE  DURABILITY 

OF  IFRA  DISTRIBUTIONS 
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TABLE  G-l.  ( CONT  *D )  SAMPLING  PLAN  FOR  M  FOR  IFRA  DISTRIBUTIONS 


■  u 


:  co 
•  uj 

CL 


Z  < 
3  U. 


D  D 

CL  OL 
UJ  UJ 
CO  CO 
£  £ 
=>  =5 
Z  Z 


ro 

v 


V 


•  Q 

z 

UJ 

*  £ 
♦  £ 

*  O 
»  u 

*  UJ 

*  a: 


132 


TABLE  G-l.  ( CONT • 0 )  SAMPLING  PLAN  FOR  N  FOR  IFRA  DISTRIBUTIONS 


■ 

•  <_> 

Z  on 

UJ 
QC 
C*  3 
H-  -» 

z  < 

3  UL 


a  a 

a;  a 

LU  1U 
tfi  CD 

SC  c 
3  3 


m 


v 


V  O 

b-« 

SC  ►- 


«  Q 

Z 


♦  a 

*  u 

*  LU 

♦  CL 


133 
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APPENDIX  H.  TABLE  H-l.  SAMPLING  PLAN  FOR  THE  MEAN 

OF  IFRA  DISTRIBUTIONS 
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TEST  TINE/DURABILITY  CRITERION 


TABLE  H-l .  CCONT*  0)  SAMPLING  PLAN  FOR  THE  MEAN  OF  IFRA  DISTRIBUTIONS 
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TEST  TIHE/DURAB ILITT  CRITERION 
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TEST  TINE/DURABILITY  CRITERION 
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TEST  TIME /DURABILITY  CRITERION 


NOTATION 


F(T) 


R(T) 

T 


T* 

a 


T 


R(a,T) 


R(«,T) 


N 

r 


C 


durability  probability  distribution  function. 
The  probability  that  a  unit  on  test  will  fail  on 
or  before  time  T  is  given  by  F(T). 

1  -  F(T)  is  the  probability  that  a  unit  on  test 
will  survive  longer  than  time  T. 

time  of  operation  of  a  unit  under  specified 
operating  conditions.  Time  is  used  in  a  general 
sense  to  represent  miles,  cycles,  number  of 
rounds,  actual  time,  or  any  other  measure  of 
operation  sequence. 

the  required  durability  life  for  an  item. 

confidence  coefficient.  If,  for  example,  a  —  .05, 
then  we  would  construct  a  100(1  -  a)%  —  95% 
lower  confidence  limit  for  R(T). 

the  actual  test  time  for  all  units  which  survive 
the  durability  test  (often  called  test  truncation 
time). 

a  100(1  -  a)%  lower  confidence  limit  for  the 
durability  life  of  an  item  relative  to  time  T. 

an  approximate  100(1  -  a)%  lower  confidence 
limit  for  the  durability  of  an  item  relative  to 
time  T. 

number  of  units  on  durability  test. 

number  of  units  which  survive  the  durability 
test. 


N-r,  the  number  of  units  which  fail  the  durabil¬ 
ity  test. 


151 


M 


f(V 

IFR 


r/T*,  is  the  normalized  test  time  for  surviving 
units,  i.e.,  each  unit  which  survives  is  on  test 
for  time  MT*.  For  example,  if  the  durability 
requirement  is  for  a  time  T*  =  5000  miles  and 
surviving  units  are  kept  on  test  for  r  =  8000 
miles,  then  M  =  8000/5000  =  1.6. 

durability  p.d.f.  =  dF(T)/d(T) 

Increasing  failure  rate 


IFRA 

R* 

r(T) 


Increasing  failure  rate  average 
Required  durability  relative  to  T* 
Failure  rate  for  a  durability  distribution 
Average  value  of  r(T) 
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